Research Article 



IJCRR 


Section: Healthcare 


PROPAGATION OF THE SHOCK WAVES 
IN THE EXPONENTIALLY DECREASING 
ATMOSPHERE 


Sci. Journal 
Impact Factor 
4.016 

Rajan Singh 1 , Mukesh Chandra 2 , Anand kumar 3 


'■*Department of Mathematics, l.F.T.M. University, Moradabad, 244001, India; 3 Ramabai Govt. Girls P.G. College Akbarpur, Ambedkamagar 
(U.P.). 


ABSTRACT 


The purpose of the present paper is to study propagation of the shock waves and the flow parameters in the exponentially 
decreasing atmosphere. Plane shock waves as it propagates vertically upward in the static atmosphere of the earth are also 
studied. The variation of Mach number and the velocity of the shock waves are obtained analytically by using Whitham’s Rule. It 
is found that shock velocity increases as the shock wave propagates in the atmosphere with decreasing density normal to plane 
of the shock front 
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1. INTRODUCTION 

Bhatnagar and Sachdev [1] have obtained a relation between 
the density, pressure and Mach number by using Whitham’s 
Rule to the propagation of the shock waves in the non- 
homogeneous medium. Kopal and Carrus [2], Hardy and 
Grover [3] have been studied the problem of propagation of 
the shock waves in a non-uniform medium with various den¬ 
sity distributions and found the behavior of the fluid flow 
in the presence of the shock waves. Hayes [4] investigated 
the vorticity jump across a gas dynamics discontinuity by 
considering the radiative heat transfer into account and con¬ 
cluded that for an optically thin gas, vorticity is unaffected in 
pseudo-stationary flows. Kanwal [5, 6] employed the theory 
of generalized functions and differential geometry to study 
the propagation and defonnation of wave front in stationary 
three dimensional gas flows. With minor changes the experi¬ 
mental data for temperature distribution given by Mitra [8,9] 
has been used .So the aim of the present paper is to see how 
the shock velocity varies in the absence of any body force 
as the shock waves propagate in the atmosphere, in which 
the density is decreasing exponentially, Whitham’s Rule [10] 
is applied to find and approximate analytic relation for the 
Mach number and the shock velocity. 

In this paper, we have studied the problem of propagation of 
the shock waves in the earth’s atmosphere and also used the 


Whitham’s rule. The experimental data for temperature dis¬ 
tribution given by Mitra [8] has been used for minor changes. 


RESEARCH METHODOLOGY 

The following Research Methodology is adopted for the pro¬ 
posed Research paper: 

• Identification of the problem 

• Collection and study of related literature 

• Mathematical formulation of the problem 

• Analysis and numerical solution of the mathematical 
model 

• Interpretation of results 

• Conclusion 

2. FORMULATION OF THE PROBLEM 

Let us assume that a gaseous medium, in which the pressure, 
the temperature and the density vary along a fixed direction. 
Let us consider that the X -axis of the medium be taken along 
this direction and the distance measured along this axis be 
denoted by x. So that the thermo dynamical parameters. In 
the direction normal to the X -axis are constant. Let the pres¬ 
sure P 0 the temperature T 0 and the density p 0 at a distance x 
measured from a fixed plane denoted by x = 0 and P , T and 
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p s be the values of the corresponding parameters at x = 0. 
Then it is considered that the density and the temperature in 
the medium are 


Po=Pse~ a ' X 

(2.1) 

T —T e ~ a 2 X 

1 o ~ 1 S t: 

(2.2) 


P = 


PO¬ 

Ps 


(2.9) 


P 


P0_ 

Ps 


( 2 . 10 ) 


And the dimensionless quantities r, R and the time t as 


Where x is the distance measured normal to the plane of vari¬ 
ations of the density and temperature. The subscript’s’ repre¬ 
sents the value of the parameters P (| ... and T at a distance x 
=0 and (X [ , CX 2 are constants. It is considered that the plane 
shock waves start at a point x = 0 and propagates in the x- 
direction. 


r = i\)(/.\ 

( 2 . 11 ) 

R = Rqcc\ 

( 2 . 12 ) 


Let us consider that the medium obeying the law of gas 


Po = R Po T o 


(2.3) 


c s { 0 a l 

" P 


(2.13) 


Where R is the gas constant and Po is the gas pressure at 
a distance x from the Surface x = 0 measured normal to it. 
Using the equations (2.1) and (2.2) into the equation (2.3), 
we obtain 


c 


s 


2 _ yp s 

Ps 


(2.14) 


The equations of motion in terms of dimensionless param¬ 
eters are given by 


P 0 = Rp s T s e ( a l +a 2) x 


Wherein 



(2.15) 


Ps = R Ps T s 


(2.5) 


dp dp du „ 

■di +u ^ +p ^ = 0 


(2.16) 


u = 


pu 0 


C 


S 


( 2 . 6 ) 


Let us consider that U 0 be the shock velocity, u 0 be the fluid 
velocity, c 0 be the sound velocity and R o be the shock radius. 
Then the non-dimensional fluid parameters are given by the 
following relations: 


p + pu 1 = P x + p x (u x —U} 2 (2.17) 

Here it is considered that there is no body force acting on the 
medium. Let a plane shock wave, created at the plane r = 0, 
propagates in the direction of r increasing. 

It is considered that the fluid velocity in front of the shock 
surface is zero, when it reaches the distance r = R. 

The jump conditions in the fluid parameters are 


4r u 0 



p 


= P0 
Ps 


(2.7) 


( 2 . 8 ) 


p + pit = P\ + /?, (zq — U) 2 

(2.18) 

p + pu 2 -Pi + p\{u i _ U) 2 

(2.19) 

v + V- »' + ‘(u,-C/) 2 

(2.20) 

r-ip 2 r- 1A 2 
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2.1 For velocity 

Substituting the equation (2.18) into the equation (2.19), we p 

obtain = y 

P 


(2.29) 


p y = p + p(Ju { (2.21) 

Further, substituting the equation (2.18) into the equation 
(2.20), we get 


= ——+ 1 ~ u2 (2.22) 

y-1 p 2 ' 

2 P 

Using the equation (2.21), the relation c = y — and U=M 
c into P 

the equation, we obtain 


{U-u,) 




u i) + 


r Pi 


y-l pu 


u\(R,t) 



y +1 


M 


(2.23) 


2.2. For pressure p,: 

1 2 ]/ 

Further: substituting U = Me and c = y — into the equa¬ 
tion (2.21), we obtain P 


M = — (2.30) 

C 

The relations (2.23), (2.25) and (2.27) are Rankine-Hugoniot 
jump conditions. 


3. SOLUTION OF THE PROBLEM: 

The aim of this chapter is to find how the shock behaves, 
when it is propagates from a plane r = 0tor=°o.For this 
it is necessary to know, how the Mach number M varies in 
the medium. This can know from the relations (2.23), (2.25) 
and (2.27), we have one extra relation connecting the four 
unknown variables u p p p p x and M in terms of p and p . 
This extra relation is provided by Whitham’s Rule [13]. The 
equation of motion along the positive characteristic axis - 

dR 

— = u, + c, 
dt 1 1 

Are given by 


Py = p + pMciiy (2.4) 

Putting the value of ul into the equation (2.24), we obtain 


P\{R,t) 

y +1 

f(M) = 

2M 2 -fl —— 
l 7) 


(2.25) 

(2.26) 


dp\ + p\c\du\ = 0 


(3.1) 


On differentiating the equations (2.23) and (2.25), we obtain 
the positive characteristic axis - 


du, 


(y+ 1) 


dp dp 

P P . 


, 1 ) 2c M 2 + 1' , 


dp - 


py 


(y+ 1 ) 


.,, 2 dM dp 

AM~ + f(M) — 

M p 


(3.3) 


2 . 3 . For density Pi : 

Inserting the equation (2.23) and U = Me into the equation 
(2.18), we obtain 


g(m) 

(2.27) 

g(M) = [{y+\)A4 2 +2\ 

(2.28) 


and differentiating the equation (2.29), we get 


2^__ dp__dp^ (34) 

C p p 

Substituting the equation (3.4) into the equation (3.2), we 
obtain 


, c (^dc \(,, 1 ) 2c f M 2 + 1) 

(V+ 1)1 cA M) (r + lA M 2 

The equation (2.29) can be written as 


dMa 


(3.5) 
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Or 


c,= P 

V A 


(3.6) 


1 


K(M)R 


Substituting the values of p t and A into the equation (3.6), 
we obtain 


M=e 2 

We shall now discuss the following two cases: 


(3.14) 


c, =■ 


-ylrf(M)g(M) 


(3.7) 


Case 1: Isothermal Case: The equation (3.13) becomes, if 
we take JJ =0 and y= 1 then 


M(y+1) 

Further, making use of the relations (2.27), (3.3), (3.5) and 
(3.7) in (3.1), we obtain 


K(M) = 


f(M) 


[2M 2 + (M 2 + 1)/?(M)] 


(3.15) 


2 [(M 2 + l)*(M) + 2M 2 ]^ + /(M) — + 2(M 2 -1)A(M)— = 0 (-3 g) 

Wherein 


If we take 77 = 1 and J = 1 then the equations (2.26), (2.28) 
and (3.9) become 

f(M) = 2M 2 ,g(M) = 2 and h(M) = M (3.16) 
Using the relation (3.16) into the equation (3.15), we obtain 


h(M)= 


g(M) 


(3.9) 


K{M) 


2 M 


(3.17) 


Hence the equation (3.8) gives the relation between M, p and 

c. 

The law of variation of the sound velocity and the pressure 
at a direction 


r = R in non-dimensional parameters is 

c(R) = e 2 

p(R) = e~ (T1+v>R 


zR 
,2 R 


(3.10) 

(3.11) 


Where 


n 


a! 

a , 


(M +1) 2 

Case II: General Case: ta) If we take 13 = 1 and y = 4/3 then 
the equation (3.13) reduces in the form 


[2M- + (. M 2 + \)h{M)\ 

we take 1 and y = 4/3 ,then the equations (2.26), (2.28) 
and (3.9) become 




Further, making use of the relation (3.19) in (3.18), we obtain 
8M 2 -1 


Inserting the relations (3. 10) and (3. 11) into the relation 
(3.8), we obtain 


K(M) = 


,„ 2 ,, 8 M -1 

- + (M- 1)J—.- 

v M 2 + 6 


2M 1 + (M 2 + 1)J - 1 

V M 2 + 6 


(3.20) 


M dR 


Wherein 


M = e 


(3.12) 


(3.13) 


(b) if we take rj= 1 and = 7/5 ,then the equations (2.26), 
(2.28) and (3.9) become 


f(M) = 2 


, ll 


M 2 ] 

[ 

M 2 -- 

7. 

>g(M) = 2 

1 + — 

5 

and h(M) = J- 


\7M--\ 
M 2 +5 


(3.21) 


Integrating the equation (3.12), we obtain 

■■e* 


^2 _ „K(M)R 


Substituting the equations (3.21) into the equation (3.18), we 
obtain 
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K(M) = 


\jM 2 — 1 ]+ (M 1 — 1)^ 


7M 1 -\ 
M 2 + 5 


2M 2 +(M 2 +1)J—^ - 1 

v M 2 +5 


(3.22) 


In the table (2.1), the function K(M) versus M is given for I] 
=0 and 7=1 and 7=1, 7=7/5 and > = 4/3. 

It is seen that the variation in K(M) as M — 2 is small lor// 
=1 and 7=1, K(M) —> 0 as M —> oo for the former case. We 
take the variation of K(M) as compared to the variations of 
M negligible. But for // = 0, K(M) varies from 0.4444 to 0 
as M goes from 2 to oo. Thus, if we take K(M) constant for 
small variation in M. The error involved is significant. The 
problem is considered in two sections, first a general case JJ 
= 1, 7 = 4/3, 7/5 and the second as isothermal case i.e. 7 

= 0,7=1- 

Table 2.1: Variation of K(M) versus M for various val¬ 
ues of 7 and 7 : 


Mach Number 

M 

^=0 

^=1 
^ = 7/5 

II 

c3 

1.0 

0.5000 

0.8570 

0.8754 

1.5 

0.4800 

1.1176 

1.1231 

2.0 

0.4444 

1.2380 

1.2365 

3.0 

0.3750 

1.3386 

1.3302 

4.0 

0.3200 

1.3775 

1.3658 

5.0 

0.2778 

1.3961 

1.3827 

6.0 

0.2449 

1.4065 

1.3920 

7.0 

0.2188 

1.4128 

1.3976 

8.0 

0.1975 

1.4169 

1.4013 

9.0 

0.1800 

1.4197 

1.4038 

10.0 

0.1653 

1.4218 

1.4056 

11.0 

0.1528 

1.4232 

1.4070 

12.0 

0.1420 

1.4244 

1.4079 

13.0 

0.1327 

1.4253 

1.4088 

14.0 

0.1244 

1.4260 

1.4094 

15.0 

0.1171 

1.4266 

1.4099 

16.0 

0.1107 

1.4270 

1.4103 

17.0 

0.1049 

1.4274 

1.4106 

18.0 

0.0997 

1.4273 

1.4109 

19.0 

0.0950 

1.4280 

1.4112 

20.0 

0.0907 

1.4283 

1.4114 

25.0 

0.0739 

1.4291 

1.4121 

30.0 

0.0624 

1.4295 

1.4124 

35.0 

0.0540 

1.4297 

1.4127 


40.0 

0.0475 

1.4299 

1.4128 

45.0 

0.0425 

1.4300 

1.4129 

50.0 

0.0384 

1.4301 

1.4130 

60.0 

0.0322 

1.4302 

1.4131 

70.0 

0.0277 

1.4303 

1.4131 

80.0 

0.0243 

1.4303 

1.4 132 

90.0 

0.0217 

1.4101 

1.4112 

100.0 

0.0196 

1.4303 

1.4132 

200.0 

0.0099 

1.4305 

1.4133 

500.0 

0.00398 

1.4305 

1.4133 

1000.0 

0.00200 

1.4305 

1.4133 

2000.0 

0.00090 

1.4305 

1.4133 

5000.0 

0.00040 

1.4305 

1.4133 

7000.0 

0.00029 

1.4305 

1.4133 

7500.0 

0.00026 

1.4305 

1.4133 

8000.0 

0.00024 

1.4305 

1.4133 



Mach number (NO 

Figure 2.1 


4. GENERAL CASE: 

Let us assume that the variation of K(M) is small as com¬ 
pared to M, then the equation (3.12) reduces in the 
form 

2 dM 

HR=- 

M K(M) 

Integrating the equation (4.1) becomes 

_ 2 |'M dM 

~ K(M) k M 
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R = - 

K(M ) 

Where M is the value of M at R=0. 

s 

Equation (4.2) can be written as 


log 


M 

M, 


K(M)R 
M = M s e 2 

Substituting the equations (3.10) and (4.3) into the equation 
(2.30), we obtain the shock velocity 


Remark 5.1: It is to be noted that 11 must be zero for isother¬ 
mal case. Hence equation (5.5) reduces in the fonn c = I and 
the equation (2.30) reduces in the form 


M=U (1.6) 

This manifests that the Mach number is equal to the shock 
velocity. By virtue of equations (5.3) and (5 (»). we obtain all 
exponential for shock velocity - 


U 2 -1 


U 2 e U 



(5.7) 


U = M s e 


n~K(M )| 


(4.4) 


Remark 4.1 : It is noteworthy that the expressions (4.3) and 
(4.4) hold only for 77 = 1, since the variation in K(M) is not 
small for 77 = 0 . 


In the fig.-2.1 and fig.-2.2, the shock velocity U has been 
drawn against the work position R for 77=1, y =\, y = 4/3, 
7/5 and 77 = 0, y =1 .It is easily seen that the acceleration 
of the shock front is small for the isothennal case compared 
to that of the general case and the strength of the shock in¬ 
creases in the medium. 


CONCLUSION 


5. ISOTHERMAL CASE: 

The variation in K(M) is not small for 1) = 0 and —>=1 
and also K(M) —> 0 as M —> 00 . It is not possible to apply 
the above method in this case. But it is easily seen that K(M) 
is very simplified for 77 = 0 , /l=l, then 


In this paper, if is very useful discussed the Plane shock 
waves as it propagates vertically upward in the static atmos¬ 
phere of the earth studied. It is found that shock velocity 
increases as the shock wave propagates in the atmosphere 
with decreasing density normal to plane of the shock front. 


K(M) = 


2 M 


[M +1 ) 2 

Using equation (5.1), then the equation (3.12) becomes 


(5.1) 


dM 

dR 


M 

M +1 




(5.2) 


Integrating (5.2), between the limit R-0 to R and M=M s to M. 
we get 
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M 2 -1 

M 2 e^ r =Ae R 

Wherein 

m}- 1 

A = M 2 s e Ms 

The expression for the dimensionless sound velocity c is 
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